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ABSTRACT 

We used the near-IR imager /spectrograph LUCIFER mounted on the Large 
Binocular Telescope (LBT) to image, with sub-arcsec seeing, the local dwarf 
starburst NGC 1569 in the JHK bands and Hel 1.08/im, [Fell] 1.64/xm and 
Br7 narrow-band filters. We obtained high-quality spatial maps of Hel 1.08//m, 
[Fell] 1.64/im and Br7 emission across the galaxy, and used them together with 
H ST / ACS images of NGC 1569 in the Ha filter to derive the two-dimensional 
spatial map of the dust extinction and surface star formation rate density. We 
show that dust extinction (as derived from the HQ;/Br7 flux ratio) is rather 
patchy and, on average, higher in the North- West (NW) portion of the galaxy 
[Eg(B-V) ~ 0.71 mag] than in the South-East [Eg(B-V) ~ 0.57 mag]. Similarly 
the surface density of star formation rate (computed from either the dereddened 
Ha or dereddened Br7 image) peaks in the NW region of NGC 1569, reaching 
a value of about 4 x 10~^ Mq yr~^ pc~^. The total star formation rate as esti- 
mated from the integrated, dereddened Ha (or, alternatively, Br7) luminosity is 
about 0.4 Mq yr~^, and the total supernova rate from the integrated, dereddened 
[Fell] 1.64/im luminosity is about 0.005 yr^^ (assuming a distance of 3.36 Mpc). 
The azimuthally averaged [FeH] 1.64/xm/Br7 flux ratio is larger at the edges of 
the central, gas-deficient cavities (encompassing the super star clusters A and B) 
and in the galaxy outskirts. If we interpret this line ratio as the ratio between 
the average past star formation (as traced by supernovae) and on-going activity 
(represented by OB stars able to ionize the interstellar medium), it would then 
indicate that star formation has been quenched within the central cavities and 
lately triggered in a ring around them. The number of ionizing hydrogen and he- 
lium photons as computed from the integrated, dereddened Ha and Hel 1.08//m 
luminosities suggests that the latest burst of star formation occurred about 4 
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Myr ago and produced new stars with a total mass of ~1.8 x 10^ Mq. 

Subject headings: galaxies: dwarf — galaxies: individual (NGC 1569) — galaxies: 
irregular — galaxies: starburst — galaxies: star formation 
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1. Introduction 

Understanding the star formation activity and history of local dwarf galaxies is of 
high astrophysical relevance. In fact, they represent ideal laboratories to study how star 
formation occurs in low-metallicity environments that, at the same time, do not show any 
global, ordered kinematic structure as in disk galaxies. As such, dwarf galaxies allow us to 
study how stellar feedback alone affects the physics, kinematics and chemical enrichment 
of the interstellar medium. Nearby, star-forming dwarf galaxies that are metal-poor and 
gas-rich are often compared to the first galaxies in the early Universe (see Izotov & Thuan 
1999), which were the first structures to collapse from primordial density fiuctuations, 
and gave rise to larger systems through mergers (White & Prenk 1991; Kauffmann, White 
& Guiderdoni 1993; Cole et al. 1994). They are also referred to as the possible local 
counterparts of star bursting galaxies that, at moderate redshift, dominate the faint galaxy 
counts at blue wavelengths (Broadhurst, EUis & Glazebrook 1992; Lilly et al. 1995). 

The nearby galaxy NGC 1569 is considered to be the archetype starbust dwarf galaxy 
for several reasons. First of all, it is gas-rich and relatively metal-poor: at the revised 
distance of 3.36 Mpc (Grocholski et al. 2008, so that = -18), its HI and dynamical 
masses are ~2 x 10* Mq and 5.1 x 10* M©, respectively (adapted from Reakes 1980). 
Its gas-phase and stellar metaUicities are estimated to be 12 -|- log(0/H) ~ 8.3 (Calzetti, 
Kinney & Storchi-Bergmann 1994; Gonzalez Delgado et al. 1997, Kobulnicky & Skillman 
1997) and Z - 0.1 - 0.2 Zq (Grocholski et al. 2008; Aloisi et al. 2001), respectively. 
Star formation has taken place in NGC 1569 for a Hubble time, although with different 
intensities as obtained by Vallenari & Bomans (1996), Greggio et al. (1998), Aloisi et al. 
(2001) and Angeretti et al. (2005) using color-magnitude diagrams of resolved stars. The 
recent studies by Angeretti et al. (2005) and Grocholski et al. (2008) show that the star 
formation activity of NGC 1569 started about 10 Gyr ago and proceeded until ~1 Gyr ago 
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presumably at a constant, very low rate. In the last 1 Gyr, however, NGC 1569 experienced 
at least three major bursts of star formation: i) the older one started about 1 Gyr ago and 
ended ~100 Myr ago with an average star formation rate (<SFR>) of about 0.04 Mq yr~^; 
ii) the intermediate episode was triggered about ~ 100 Myr ago and lasted up to 30 Myr 
ago with <SFR> ~ 0.08 M© yr~^, and in) the younger burst begun ~ 27 Myr ago and 
ended about 8 Myr ago with <SFR> ~ 0.3 Mq yr~^ (Angeretti et al. 2005). According to 
Aloisi et al. (2001), the spatial distribution of stars across the galaxy changes with their 
age, so that stars younger than 50 Myr are more centrally concentrated, intermediate age 
stars (50 Myr - 1 Gyr) are uniformly distributed, and older stars are mostly found in the 
outskirts of NGC 1569. Angeretti et al. (2005) suggested that this episodic star formation 
(in the last 1 Gyr) could be ascribed to the gravitational interactions of NGC 1569 with an 
HI cloud (~ 7 X 10^ Mq in mass) a few kpc away yet connected to the galaxy by an HI 
bridge (cf. Stil & Israel 1998). 

The presence of several HII regions (Waller 1991) and the spatially extended Ha 
emission indicate that NGC 1569 has formed new stars rather recently. Hunter & Elmegreen 
(2004) used the spatially integrated Ha luminosity (corrected for reddening) to derive 
a SFR of ~ 0.6 M© yr~^ (adjusted to a distance of 3.36 Mpc). The morphology of the 
extended Ha is rather complex and characterized by arcs, filaments and four large-scale 
superbubbles (Hunter et al. 1993; Heckman et al. 1995; Martin 1998; Westmoquette et 
al. 2008). The latter are seen to expand at typical velocities of ~100 km s~"^ which imply 
dynamical ages < 25 Myr, consistent with the most recent burst of star formation derived 
by Angeretti et al. (2005; Heckman et al. 1995; Westmoquette, Smith & Gallagher 2008). 
The expansion of these superbubbles is most likely powering a galactic outflow throughout 
the disk whose direction is approximately perpendicular to the inclined and flattened HI 
disc of NGC 1569. The presence of an outflow of ionized gas is also supported by the high 
[SII]/Ha ratios measured within the superbubbles and indicative of shocks (Westmoquette, 
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Smith & Gallagher 2008). The arcs of ionized gas detected in Ha and X-ray in the 
super bubbles suggest that the hot gas is still confined within the super bubbles, although it 
is moving at the escape velocity of NGC 1569. Thus, at present day, the galactic outflow of 
NGC 1569 does not appear to be in a steady state and to be able to induce mass loss from 
the galaxy (Westmoquette, Smith & Gallagher 2008). 

NGC 1569 is special also for another reason: it hosts two of the nearest super-star 
clusters known (SSCs A and B, Arp & Sandage 1985). SSCs are massive star clusters 
few Myr to several hundred Myr old, believed to be the ancestors of globular clusters. In 
fact, they are as extended and massive as globular clusters and, after ~10 Gyr of passive 
evolution, they can match the luminosity range spanned by globular clusters today (cf. van 
den Bergh 1995; Meurer 1995). They are typically observed in interacting/merger systems 
(see Whitmore et al. 1993) and starburst galaxies (see Meurer et al. 1995). In NGC 1569, 
SSC A has been resolved into two components (Al and A2, De Marchi et al. 1997), of 
which A2 may be the host of young O and Wolf-Rayet stars and Al together with SSC B 
may be dominated by older red supergiants (Origlia et al. 2001). Gonzalez Delgado et al. 
(1997) suggested that SSCs A and B have possibly undergone sequential star formation, 
with two major bursts of star formation about 3 and 9 Myr ago. The authors also found 
a deficit of ionized gas around SSCs A and B which may have been created by the stellar 
winds and supernova explosions of the older burst removing the local gas. According to the 
measurements of Ho & Filippenko (1996) SSC A is as massive as ~4 x 10^ Mq (adjusted 
to a distance of 3.36 Mpc). 

Age dating of star clusters and resolved stars strongly depends on dust extinction. 
Studies of the stellar content of NGC 1569 based on long-slit spectroscopy show that the 
intrisic dust extinction varies across NGC 1569 by few tenths of a dcx, and peaks at SSCs A 
and B (Gonzalez Delgado et al. 1997; Origlia et al. 2001) and in the North- West portion of 
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the galaxy (Kobulnicky & Skillman 1997). The spatial coverage of long-slit spectroscopy is 
coarse, however, and no uniform and contiguous two-dimensional map of the dust extinction 
in NGC 1569 is available in the literature. Clearly, such a 2D map helps in reaching a 
higher accuracy in the determination of the star formation history and its associated star 
formation rate across NGC 1569. For this reason, we have taken advantage of the large 
field of view of LUCIFER (LBT NIR spectroscopic Utility with Camera and Integral-Field 
Unit for Extragalactic Research) mounted on the Large Binocular Telescope (LBT) to 
perform deep imaging in the Br7 filter and in the K band. Our aim is to construct the 2D 
spatial map of dust extinction across NGC 1569 from the Br7 and Ha (obtained with the 
Advanced Camera for Surveys on HST) images of the galaxy and to derive the 2D spatial 
distribution of star formation rate density. We also observed the Hel 1.08// and [Fell] 1.64// 
lines and their broad-band continua in order to estimate the number of massive young stars 
and the supernova rate in NGC 1569, respectively. The observations and data reduction are 
described in Sect. 2. The spatial distribution of the color excess Eg(B-V) associated with 
the gas and the spatial distribution of the star formation rate density across NGC 1569 are 
derived in Sect. 3 and 4, respectively, while the properties of the [Fell] 1.64/fm emission are 
analyzed in Sect. 5. Conclusions follow in Sect. 6. 



2. Observations 

The observations presented here were carried out with the LBT, located on Mount 
Graham, Arizona (Hill, Green & Slagle 2006). NGC 1569 was observed during the science 
commissioning of LUCIFER (Ageorges et al. 2010; Seifert et al. 2010), with the N3.75 
camera which provides a 4'x 4' field of view at an angular resolution of 0.12"/pixel 
(equivalent to 1.95 pc per pixel at the adopted distance of 3.36 Mpc). We imaged NGC 1569 
between September and November 2009 in the JHK broad bands and in the narrow-band 
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Hel 1.08/im, [Fell] 1.64/im and Br7 filters, with a seeing < 0.5" (i.e. 8 pc) throughout 
the observing runs. Although the exposure times in the broad-band and Br7 filters were 
estimated from the available 2MASS photometry and H ST /ACS F658N (Ha) image via 
the LUCIFER Exposure Time Calculator, the effective exposure times were largely dictated 
by the available amount of time during the LUCIFER commissioning. The exposure time 
for each of the Hel 1.08/im and [Fell] 1.64/im filters was initially set to 1 hour (following 
Labrie & Pritchet 2006), but it was later cut down to 30 min for the Hel 1.08/im. The final 
exposure times on source achieved for the near-IR dataset are listed in Table 1, together 
with the imaging taken in the F606W (broad V) and F658N (Ha) filters with HST/ ACS 
as part of program ID 10885 (PI Aloisi). 

The near-IR images were reduced with standard IRAlj^ routines. They were corrected 
for bias and fiat-fielding, and their astrometric solution was derived from the field stars in 
common with 2MASS. After background subtraction, the images taken in the same filter 
were corrected for geometric distortion and then combined together on the basis of their 
World Coordinate System with a weighted mean. The final mosaiced frames in the broad 
band filters were fiux-calibrated, at a 10% accuracy, using the 2MASS magnitudes of the 
stars in the field of NGC 1569. 

After convolution to the same PSF (FWHM = 0.5") the continuum emission was subtracted 
from each narrow-band image by scaling the closest broad-band image (J for Hel 1.08/im, 
H for [Fell] 1.64/im and K for Br7) so that field stars would have had the same integrated 
count rate in both the narrow and broad band filters. This scaling was applied to the broad 
band zero points in order to calibrate the pure emission line images in fiux, with an achieved 

^IRAF is distributed by the National Optical Astronomy Observatory, which is oper- 
ated by the Association of Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 
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accuracy of 20% (estimated from the scatter in the zero point derived from different stars). 

The optical HST/ACS images were retrieved from the HST archive already reduced. 
The subtraction of the continuum emission (the F606W image) from the F658N frame 
and the flux calibration of the pure Ha emission were performed as for the near-IR data, 
reaching a flux accuracy of 20%. Such a low accuracy in the flux calibration of the Ha is 
mostly due to the broad F606W filter which makes it impossible to account for color terms. 
The Ha-emission image was corrected for a 5% fiux contribution from the [Nil] A6548,6584 
lines (Moustakas & Kennicutt 2006), and then convolved with a Gaussian function to have 
the same stellar PSF as the Br7-emission image, resampled to achieve the same pixel scale 
as the LUCIFER N3.75 camera, and aligned to the Br7-emission image. 

Most of the analysis presented in the following sections requires or involves the 
production of spatially resolved maps of fiux ratios in different passbands. In order to 
ensure a sufficient S/N for the largest number of pixels while retaining the best possible 
spatial resolution, we adopt the approach introduced by Zibetti, Chariot & Rix (2009) using 
the adaptive smoothing code ADAPTSMOOTH (Zibetti 2009). We first smooth each image 
individually by taking the median fiux among the pixels inside a top-hat circular kernel, 
whose size is determined as the minimum that allows to obtain a S/N> 10 (background 
dominated noise is assumed). The spatial match between different bands is obtained by 
requiring that the kernel size at each pixel location ensures a minimum S/N of 10 in all 
images (see Zibetti 2009 for details). The smoothing procedure also results in a first S/N 
cut: all pixels for which sufficient S/N of 10 cannot be reached in all images even when 
the maximum kernel size of 27 pixels is used are fiagged. A further S/N cut is apphed by 
excluding all pixels whose fiux (in the smoothed image) is less than 3a of the large scale 
background fluctuations as measured from randomly placed boxes of a few hundred pixels 
on a side. 
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3. The spatial map of dust extinction 

The hne-emission images of NGC 1569 are shown in Fig. 1, aheady flux-cahbrated 
but not corrected for dust extinction. Fluxes, in units of erg cm~^ s~^, are displayed in 
logarithmic scale across an area of 1.5 kpc x 1 kpc at the adopted distance of 3.36 Mpc. We 
clearly detect the base of the north-western arm as well as "cavities" at the center of the 
galaxy where SSCs A and B reside, and corresponding to the cluster surroundings where 
Gonzalez Delgado et al. (1997) found a deficit of ionized gas. These regions were flagged as 
underflow pixels by our smoothing procedure. The S/N cut used in our image smoothing 
results in a hmiting flux of 10"^^, 10~^^-^, pa 10~^^ and pa 10"^^-^ erg cm""^ s~^ for Ha, 
Hel 1.08/im, [Fell] 1.64/im and Br7, respectively (corresponding to the lower limits of the 
flux ranges displayed in Fig. 1). The morphologies of the Ha and Hel 1.08/xm emissions, 
which trace the location of young and massive stars, are rather similar. They exhibit two 
major peaks along the major axis of NGC 1569, one on each side of the cavities and with 
the brighter peak located in the North- West (NW) half of the galaxy. Relatively high 
fluxes are also detected along nearly the full perimeter of the cavities. The Br7 emission is 
morphologically consistent with that of the Ha although it seems to be relatively weaker 
along the eastern rim of the cavities. Finally, the [Fell] 1.64/x emission, which is typically 
due to shocks induced by supernova explosions (but see Sect. 5), is as diffuse as the H and 
He emissions. The [Fell] flux peaks only in the NW half of NGC 1569 and stays relatively 
high along the southern rim of the cavities. The very bright compact source in the eastern 
part of NGC 1569 is a supernova remnant (cf. Labrie & Pritchet 2006). 

We integrate the observed Ha and Br7 maps of NGC 1569 over an eUipse of semi-major 
axis a = 37" (606 pc) and semi-minor axis b = 17" (274 pc) which does not include the 
extended Ha superbubbles, where shocks due to the galactic outflow seem to be important. 
The depth of our near-IR imaging is not sufficient to allow us to detect these extended 
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superbubbles. Pixels below our S/N cut are excluded from the integration^ We use 
the resulting total Ha and Br7 fluxes to derive the color excess (intrinsic + foreground) 
averaged across the galaxy as in Calzetti et al. (1996): 

^ -R m ^n^_ -log(Robs/Rint) .y. 

< -tvglO - Vj >- 7 -T TT (Ij 

0.4[K(AHa) - /«(ABr7)J 

where Robs is the observed Ha/Br7 flux ratio, Hint is the intrisic Ha/Br7 flux ratio (102.8) 
computed for Case B (i.e. the gaseous clouds are optically thick to line photons) assuming 
a temperature of 10,000 K and a density of 100 cm~^ (Osterbrock 1989). The extinction 
curve is adopted from Fitzpatrick (1999), with /€(Ahq) = 2.535 and k^Xbt-y) ~ 0.363. 
We thus obtain < Eg(B — V) > = 0.58 mag, consistent with Calzetti et al. (1996), and an 
intrisic < Eg(B — V) > = 0.07 mag assuming a foreground Galactic color excess Eg(B — V) 
= 0.51 mag (Burstein & Heiles 1982; Calzetti et al. 1996). Since the intrinsic color excess 
on the stellar continuum is < Eg(B — V) > = 0.44< Eg(B — V) >, we estimate the total 
(intrinsic + foreground) color excess on the stellar continuum to be < Es(B — V) > = 0.54 
mag, consistent with the values derived by Angeretti et al. (2005) and Grocholski et al. 
(2008) from the color-magnitude diagram of resolved stars. 

We apply Eq. (1) to each pixel in the observed Ha and Br7 maps in order to derive 
the spatial distribution of Eg(B-V) (intrinsic + foreground) across NGC 1569. This is 
shown in Fig. 2; we clearly see that the color excess is higher (> 0.6 mag) where the H and 
He emission fluxes are higher, and the NW half of the galaxy is more extinguished than 
the South-East [where Eg(B-V) < 0.6 mag]. The distribution of the pixel color excess as 
derived from Fig. 2 peaks at about Eg(B-V) = 0.6 mag and its full width half maximum 



^We compared the resulting total fluxes with those obtained from the same area prior 
applying ADAPTSMOOTH to flnd negligible differences, at the 1% - 2% level. 
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points to variations in Eg(B-V) up to 0.4 mag across the galaxy. The presence of a gradient 
in Eg(B-V) along the galaxy major axis is evident in Fig. 3, where we plot the azimuthal 
Eg(B-V) as a function of position along the semi-major axis (solid line; the dotted line 
represents the foreground Eg(-B — y) = 0.51 mag). We integrate the observed flux of the 
Ha and Br7 images over concentric, elliptical semi annuli and compute the HQ;/Br7 flux 
ratio to estimate the azimuthal Eg(B-V) as a function of position along the semi- major axis 
of NGC 1569. From now on, we assume the convention of negative distances for increasing 
right ascension from the galaxy center. The vertical errorbar (±0.1 dex) in Fig. 3 indicates 
the systematic uncertainty on Eg(B-V) due to the uncertainty on the flux calibration. 
Because of our S/N cut, the color excess can be determined only for galactocentric distances 
<-50 pc and >30 pc. As Fig. 3 shows, the color excess increases from Eg(B-V) ~ 0.57 
mag in the South-East portion of the galaxy to Eg(B-V) = 0.7 - 0.8 mag in the NW. 
This gradient is consistent with the findings of Kobulnicky & Skillman (1997). The higher 
Eg(B-V) (~0.7) seen at -1-100 pc from the central cavities in Fig. 3 is in agreement with the 
findings of Gonzalez et al. (1997). 

We use the 2D spatial map of Eg(B-V) to correct the hne-emission images pixel by 
pixel following the formula: 

Fi{\) = F„(A) X ioO-4Eg(B-v)«(A) (2) 

where Fi{X) and Fo{X) are the intrinsic and observed line flux, respectively, and k{X) is 
the extinction law. For a Fitzpatrick's (1999) extinction curve k is 2.535, 1.096, 0.562 
and 0.363 for Ho;, Hel 1.083/xm, [Fell] 1.64/xm and Br7, respectively (Calzetti, private 
communication). The correction for dust extinction enables us to correct the [Fell] 1.64//m 
emission for the contribution of the Brl2 line, estimated to be 0.17 times the Br7 flux 
(assuming Case B, Storey & Hummer 1995). We then compute the surface brightness 
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profile of NGC 1569 for concentric, elliptical, semi annuli of fixed ellipticity (b/a = 0.45) 
and fixed orientation of the semi- major axis (P. A. = 117°) in each filter. These profiles 
are shown in Fig. 4 where the surface brightness is corrected for reddening but not for the 
galaxy inclination. The gap in the line-emission profiles between -50 and 30 pc is due to our 
S/N cut masking out the central gas-deficient regions. The surface- brightness profiles of the 
H and He lines share a similar shape, with maxima at a galactocentric distance of about 
±300 pc. At smaller galactocentric distances the surface brightness of the H and He lines 
declines, while that of the [Fell] line stays constant and is ~0.6 dex brighter than Br7. At 
galactocentric distances <-300 pc and >300 pc all five profiles are close to an exponential 
disk with similar scale lengths. All profiles appear to be steeper in the NW half of the 
galaxy with a typical, exponential scale length of 70 pc against 90 pc as measured for the 
SE half. The peak at -550 pc in the [Fell] surface brightness is due to the emission from a 
compact supernova remnant. 

4. The recent star formation activity of NGC 1569 

We integrate the dereddened Ha, Hel 1.08/xm and Br7 maps of NGC 1569 over an 
ellipse of semi-major axis a — 37" (606 pc) and semi-minor axis b — 17" (274 pc) which does 
not include the extended Ha superbubbles. Pixels below our S/N cut are excluded from 
the integration. We obtain a total, dereddened luminosity of 5.31 x 10"^*^ erg s^^ (±0.20 
dex), 4.51 X 10^^ erg s"^ (±0.14 dex) and 5.05 x 10^^ erg s"^ (±0.11 dex) for the Ha, 
Hel 1.08//m and Br7 lines respectively, at the adopted distance of 3.36 Mpc. The quoted 
uncertainties (all at la level) are due to flux calibration and reddening correction. We use 
these values to study the stellar population responsible for ionizing the interstellar H and 
He, and the dereddened Ha and Br7 images to trace the spatial distribution of the SFR 
density across NGC 1569. 
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4.1. The recent stEirburst 

We compute the relation between the number of ionizing photons [(5(H°)] for hydrogen 
and the Ha (617) hne luminosity, as well as the relation between the number of ionizing 
photons [(5(He°)] for helium and the Hel 1.08/xm line luminosity, under Case B with a 
temperature of 10,000 K and a density of 100 cm~^ (typical of HIT regions): 

g(H°) [s-^] = 7.3 X 10^^ L(Hq;) [erg s"^] (3) 
g(H°) [s-^] = 7.5 X 10^^ L{Bt^) [erg s"^] (4) 

g(He°) [s-^] = 1.0 X 10^2 L(HeI 1.08/xm) [erg s"^] (5) 

Applying Eq.(3) - (5) to the total, dereddened luminosities from above yields Q(H°) = 
3.88 X 10^2 s"^ from the Ha line, Q(H°) = 3.79 x 10^^ s"^ from the Br7 line, and g(He°) 
= 4.51 X 10^^ s~^. In Fig. 5 we compare Q(H°) (as derived from the Br7 line) and (5(He°) 
with STARBURST99 (Leitherer et al. 1999) predictions for an instantaneous burst of star 
formation (panel a) and for continuous star formation (panel b). For both scenarios we 
adopted Z — 0.2 Zq and a Kroupa (2001) initial mass function (with exponents = 1.3 and 
2.3 for the mass ranges 0.1 - 0.5 Mq and 0.5 - 100 respectively); stellar evolution was 
modeled with the Geneva evolutionary tracks for high mass- loss rates (Meynct et al. 1994). 
In panel a of Fig. 5 the correlation between (5(H°) and (5(He°) for different stellar ages and 
at fixed total mass (Mtot) of the burst is traced by solid lines corresponding to Mtot = 10^, 
10^ and 10^ M©. The dashed lines represent the correlation between Q{li^) and (5(Hc'^) 
for different Mtot and at fixed stellar age (1, 3, 5, and 7 Myr). The number of H and He 
ionizing photons measured for NGC 1569 from its Br7 and Hel 1.08//m lines is plotted with 
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a black-filled circle whose errorbars (±l(j) take into account the systematic uncertainty due 
to flux calibration and reddening correction. This number is consistent with being produced 
by a 4 Myr (±0.3 Myr at la level) old stellar population, comprising ~5400 O stars and 
~250 Wolf-Rayet stars and with a total mass of 1.8 x 10^ M© (with a 20% uncertainty, 
at 1(7 level). In panel b the observed values of Q(H°) and (5(He°) are compared with the 
predictions for continuous star formation: the solid lines show the correlation between 
Q(H°) and g(He°) at fixed SFR of 0.05, 0.10, 0.15 and 0.20 M© yr-^ and their length 
spans the time interval between 1 Myr and 1 Gyr. At stellar ages older than 1 Gyr, (5(H°) 
and Q(He°) cease to change with time. None of the tracks fits within la the number of 
H and He ionizing photons measured for NGC 1569; the observed (5(H°) and (5(He°) are 
marginally consistent with a star formation history with SFR > 0.15 M© yr~^ that has 
been going on for more than 10 Myr. 

We checked whether the properties of the stellar population ionizing the H and 
He gas in NGC 1569 vary with galactocentric distance. We integrated the fiux of the 
dereddened Ha and Hel 1.083;um images over concentric, elliptical semi annuli and derived 
the corresponding (5(H°) and (5(He°) as a function of galactocentric distance. We compared 
these (5(H°) and (5(He°) values with the SB99 predictions for instantaneous bursts, and 
found that stellar ages are consistent with a value of ~4 Myr at any galactocentric distance. 
On the other hand, the mass surface density of the ionizing stars increases inward, from 
~0.01 Mq pc-2 at ±700 pc from the galaxy center to ~0.3 M© pc-^ at 300 pc in the SE 
portion of NGC 1569 and to ~1.8 M© pc^^ in the North- West. 

4.2. The stair formation rate and its density 

Using the relation between star formation rate (SFR) and the dereddened Ha (or 
Br7) luminosity as in Kennicutt (1998) we derive an average SFR of ~0.4 Mq yr"^. 
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consistent with what can be inferred from the age and total mass of the burst derived 
above. The uncertainty on the flux cahbration of both emission hnes and on the reddening 
determination gives rise to a systematic uncertainty of ±0.20 dex on the SFR derived from 
the HcK (or 0.11 dex on the SFR estimated from the Br7). 

The 2D map of the star formation rate per pixel, obtained from the dereddened Br7 
image, is shown in Fig. 6 in logarithmic units; the SFR across NGC 1569 is typically ~10~^ 
M0 yr~^ per pixel (i.e. 3 x 10"^ M0 yr~^ pc~^) and increases to few 10~^ M0 yr~^ per 
pixel (i.e. ~ 10^^ Mq yr~^ pc^^) at each side of the cavities where the H and He line 
emissions peak. This trend is well described by the azimuthal SFR density (obtained by 
integrating the dereddened flux of the Br7 image over concentric, elliptical semi annuli) 
plotted as a function of position along the semi-major axis (see Fig. 7). The maxima of the 
SFR density are about 9 x 10"'^ M© yr"^ pc"^ and 4 x 10"^ M© yr"^ pc"^ at -300 pc and 
-1-300 pc from the galaxy center, respectively. The errorbars (±l(j) in Fig. 7 indicates a 
systematic uncertainty of ±0.1 dex on the SFR derived from the dereddened Br7 image. 

The width of the distribution of the pixel star formation rate as obtained from the 
dereddened Br7 image indicates that SFR density varies by a factor of 10 within the galaxy. 

5. The [Fell] emission 

The [Fell] emission, with its prominent lines at 1.26 //m and 1.64 fim, is usually 
explained as triggered by electron collisions occurring in a zone of partially ionized hydrogen 
where Fe"*" and e~ coexist. The emission intensity is proportional to the size of this region, 
which is rather small in HIT regions but extended in supernova remnants (SNRs, cf. Oliva, 
Moorwood & Danziger 1989). The [Fell] emission would thus be a clear signpost for SNRs, 
and as such it is expected to come from compact sources in external galaxies. Recent 
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studies have shown that the [Fell] emission can also be spatially extended in galaxies known 
to experience a galactic wind (e.g. NGC 253, Forbes et al. 1993, Sugai, Davies & Ward 
2003; NGC 5253, Cresci et al. 2010, Labrie & Pritchet 2006, Turner, Beck & Ho 2000; 
M 82, Greenhouse et al. 1997, Heckman et al. 1987). Some theoretical computations (e.g. 
Scab & ShuU 1983 and McKee et al. 1987) have shown that this extended [Fell] emission 
can be explained by high-speed shocks such as those associated with a galactic outflow. 
The proposed mechanism assumes that iron atoms can be extracted from silicate grains 
that are preferentially destroyed via nonthermal sputtering and grain-grain collisions (cf. 
also McCarthy, Heckman & van Breugel 1987; van der Werf et al. 1993. We note that 
these mechanisms can also be triggered by SN shocks). NGC 1569 is yet another starburst 
galaxy exhibiting spatially extended emission of [Fell] 1.64//m (already detected by Labrie 
& Pritchet) and a galactic outflow (Heckman et al. 1995). 

Given that the Br7 emission comes from gas ionized by OB stars, the flux ratio [Fell] 
1.64/im/Br7 is commonly used to trace the number ratio of supernovae (SNs) to OB stars, 
hence the star formation history of a galaxy at fixed initial mass function. Values of this 
ratio larger than few tens indicate that excitation is mostly produced by SN shocks, while 
values << 1, typically 0.1 or lower, are due to photoionization (Alonso-Herrero et al. 1997). 
The dereddened [Fell] 1.64/im luminosity (already corrected for the emission of the Brl2 
line) integrated over an ellipse of semi-major axis a = 37" and semi-minor axis 6 = 17" 
(excluding the central cavities) is 4.67 x 10^^ erg s"^ at the adopted distance of 3.36 Mpc. 
From this and the integrated luminosity of Br7, we obtain a global [Fell] 1.64//m/Br7 ratio 
= 0.9, suggesting that excitation by SN shocks may prevail over photoionization by OB 
stars formed during the most recent burst. Figure 8 shows the two-dimensional spatial 
distribution of the [Fell] 1.64//m/Br7 ratio across NGC 1569. The line ratio is larger than 
1 mainly in the galaxy outskirts and around the central cavities, and it is < 1 especially in 
two regions at the opposite extremes of the cavities coincident with the peaks of the Ha 
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emission and the star formation rate density. Another way of looking at the spatial variation 
of the [Fell] 1.64/im/Br7 ratio is to integrate the dereddened flux of the [Fell] 1.64//m and 
Br7 images over concentric, elliptical semi annuli and compute the [Fell] 1.64/im/Br7 flux 
ratio as a function of position along the semi-major axis. This is plotted in Fig. 9; the 
minima at [Fell] 1.64//m/Br7 < 1 occur at about ±300 pc from the galaxy center where the 
star formation rate density is highest (cf. Fig. 8). The peak at -550 pc is due to a compact 
SNR. The [Fell] 1.64//m/Br7 ratio increases above 1 with galactocentric distance and at 
the edges of the cavities, possibly suggesting that in these regions star formation has been 
somehow quenched. A relevant caveat on this last interpretation comes from the possibility 
that iron has been excited by shocks produced by the galactic outflow of NGC 1569. In this 
case, the [Fell] 1.64//m/Br7 ratio may be expected not to trace the star formation history 
of a galaxy on small spatial scales. The extent by which a galactic outflow can alter this 
ratio so that it is no longer representative of the local number ratio SNs vs OB stars is 
still to be thoroughly investigated. In their spectroscopic study of the central 20" x 20" in 
NGC 1569 Westmoquette et al. (2007) found little evidence for shocked line ratios, most 
likely because this area is experiencing intense star formation and photoionization from OB 
stars. In this region (equivalent to about 320 x 320 pc^ at the adopted distance of 3.36 
Mpc) we measure the lowest values of [Fell] 1.64//m/Br7 (see Fig. 9, exception made for 
the cavities edge), which may also be interpreted as due to photoionization. 

We used the dereddened, integrated luminosity of the [Fell] 1.64//m and Br7 hues to 
estimate the supernova rate (SNR) and the gas-phase abundance of Fe+ for NGC 1569 as 
a whole. Under the assumption that the typical [Fell] 1.64^m luminosity of a supernova 
is ~10^^ erg s~^ and its average life ^10"^ yr (cf. Lumsden & Puxley 1995), we estimate a 
SNR of ~0.005 yr~^ for the whole galaxy. As for the gas-phase abundance of Fe"*", we adopt 
Greenhouse et al.'s (1997) prescription, where: 
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Fe+ FeII l.64 

under Case B for a temperature of 10,000 K and a density of 100 cm~^. For the same 
physical conditions, L(Pa/3) = 5.86 L{Byj), hence: 

N(Fe+) fiL([FeII11.64) 

We thus derive N{Fe+)/N{E.) > 1.4 x 10"^ and A^(Fe)/A^(Fe)0 > 0.05, where the Sun 
A'"(Fe)/A'"(H) is 2.8 x 10~^ (Holweger 2001). The lower hmit comes from the assumption 
that all the iron is singly ionized. 



6. Conclusions 

The wide field of view and high angular resolution of the infrared imager / spectrograph 
LUCIFER mounted on the Large Binocular Telescope allowed us to acquire deep and 
detailed images of the local starburst NGC 1569 in the Hel 1.08//m, [Fell] 1.64/im and Br7 
light. Together with H ST /ACS Ha images of the galaxy, these data were used to derive 
the two-dimensional spatial distributions (on scales as small as 2 pc) of dust extinction 
(foreground + intrinsic), star formation rate density and [Fell] 1.64//m/Br7 ratio, and 
to estimate the age and total mass of the most recent burst responsible for ionizing the 
interstellar medium (ISM) in NGC 1569. The galaxy looks clearly asymmetric along its 
major axis with respect to its center, where its super-star clusters A and B have most likely 
evacuated gas and formed gas-deficient cavities. The color excess Eg(B-V), as derived from 
the HQ;/Br7 ratio, is patchy with a scatter of about 0.4 mag across the galaxy. Reddening 
is lower in the SE portion of NGC 1569 with an average Eg(B-V) ~ 0.57 mag and increases 
to 0.7 - 0.8 mag in the NW along the galaxy major axis. Similarly to the surface brightness 
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profiles of the H and He line emissions, the star formation rate density increases inward 
along the major axis, reaching two maxima of 9 x 10~^ and 4 x 10~^ Mq yr^^ pc~^ at -300 
and +300 pc, respectively, from the galaxy center (i.e. just outside the central cavities). 
Within the galaxy the star formation rate density is seen to vary by a factor of 10. The 
peaks in SFR density are spatially coincident with the minima in the [Fell] 1.64//m/Br7 
flux ratio, suggesting that the OB stars formed during the most recent burst are responsible 
for the ISM photoionization and are preferentially located in a ring around the central 
cavities. The [Fell] 1.64//m/Br7 ratio is seen increasing outward from the galaxy center and 
along its major axis, as well as at the edges of the central cavities. At face value, this trend 
may imply that excitation by SN shocks dominates over photoionization from OB stars and 
indicate that star formation has been quenched some Myrs ago in the central cavities and 
in the outermost regions of the galaxy. This trend would be consistent with the findings of 
Aloisi et al. (2001), whereby the age of stars increases with galactocentric distance from 
50 Myr close to the galaxy center to >1 Gyr in the galaxy outskirts. How and to which 
extent the intervention of a galactic outflow as in NGC 1569 alters the [Fell] 1.64//m/Br7 
ratio and prevents it from tracing the ISM excitation mechanisms on small spatial scales 
remains to be fully investigated. On global scales, the [Fell] 1.64/xm luminosity integrated 
over the galaxy main body points to a SN rate of about 0.005 yr~^ and an iron abundance 
A'"(Fe"^)/A'"(H) > 1.4 x 10"^. This value is quite consistent with the Fe abundance in the 
ISM of the Small Magellanic Cloud (Rollcston ct al. 2003). If the high-speed shocks in 
the galactic outflow would be responsible for extracting iron atoms from silicate grains (as 
suggested by the theoretical calculations of, for example. Scab & ShuU 1983 and McKee 
et al. 1987) and hence for a spatially extended [Fell] emission, one would expect to find 
the ISM Fe abundance in NGC 1569 to be enhanced with respect to the Small Magellanic 
Cloud (SMC) which has a stellar metaUicity similar to NGC 1569 but no galactic outflow. 
Given that the ISM Fe abundance is similar in both NGC 1569 and the SMC, the present 
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data would tend to exclude that the diffuse [Fell] emission in NGC 1569 is due to the 
galactic outflow. Moreover, Sofia et al. (2006) found that iron atoms in the ISM of the 
SMC are most likely locked in metal grains and/or oxides rather than in silicates. If this 
apphed to NGC 1569 also, the galactic outflow would not be able to extract iron from dust 
grains (as in the scenario proposed by Scab & ShuU 1983 and McKee et al. 1987), increase 
the Fe emission and thus explain its spatial extension. 

We used the Ha, Br7 and Hel 1.08//m luminosities integrated over NGC 1569 to 
estimate the star formation rate of the most recent burst (~0.4 Mq yr~^) and the number 
of stellar photons able to ionize the interstellar H and He gas [i.e. (5(H°) and (5(He°)]. 
When compared to STARBURST99 predictions the observed (5(H°) and (5(He°) enable us 
to constrain the age and the mass of the stellar population responsible for ionizing H and 
He. In the case of NGC 1569 the values of (5(H°) and (5(He°) are in good agreement with 
an instantaneous burst of star formation. A scenario where star formation is continuous and 
takes place with a constant SFR is not rejected by statistics, but is seemingly inconsistent 
with the gasping nature of the star formation activity undergone by NGC 1569 in the last 
1 Gyr, as derived by Angeretti et al. (2005) from the color-magnitude diagram of resolved 
stars. The best-fitting instantaneous burst occurred about 4 Myr ago and involved a total 
stellar mass of ~1.8 x 10^ Mq. It gave birth to about 5400 O stars and 250 Wolf-Rayet 
stars. Such a stellar age is constant with galacto centric distance and where Ha, Br7 and Hel 
1.08/xm emissions are detected. The only property of the burst changing with distance from 
the galaxy center is the surface mass density of the ionizing stars, which follows the surface 
brightness profiles of the H and He lines and increases from ~0.01 Mq pc^ in the galaxy 
outskirts up to ~1.8 Mq pc^ just outside the central cavities. The comparison of this latest 
starburst with the three episodes of star formation identified by Angeretti et al. (2005, cf. 
Sect. 1) suggests a star formation history where the burst duration and average SFR have 
become shorter and larger, respectively, from 1 Gyr ago to 4 Myr ago. Consequently, the 
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total mass in stars produced by these bursts of star formation has decreased by a factor of 
about 20 in the last 1 Gyr. Interestingly enough, the 4 Myr age derived from the observed 
(5(H°) and (5(He°) is close to that of the younger stellar component found mostly in SSC 
A by Gonzalez Delgado et al. (1997). Therefore, the big picture emerging from these 
results is one where the strong stellar winds and supernova explosions from the older stellar 
population of SSCs A and B (about 9 Myr old) removed a large fraction of gas from the 
clusters surroundings and triggered, ~4 Myr later, star formation at the edges of the central 
cavities. The gas left over in the vicinity of SSC A (and perhaps also B) underwent star 
formation nearly at the same time. 

We would hke to thank Daniela Calzetti and Monica Tosi for very useful discussions. 
The Dark Cosmology Centre is funded by the DNRF. 

Facilities: LBT (LUCIFER), HST (ACS) 
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Fig. 1. — The line-emission images of NGC 1569 calibrated in flux in units of erg cm~^ s^^, 
but not corrected for reddening. Fluxes are on logarithmic scale, and the image size is 1.5 
kpc X 1 kpc at the adopted distance of 3.36 Mpc. The white frames indicate the location of 
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Fig. 2. — The Eg(B-V) (intrinsic + foreground) map of NGC 1569 derived from the HQ!/Br7 
flux ratio. The image size is 1.5 kpc x 1 kpc. 
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Fig. 3. — The azimuthal Eg(B-V) (intrinsic + foreground), computed over 180 degrees, as 
a function of position along the semi-major axis (sohd hnc). The dotted hne represents the 
foreground Galactic color excess in the direction of NGC 1569, Eg(B-V) — 0.51 mag. The 
errorbar corresponds to a systematic uncertainty on Eg(B-V) of ±0.1 mag, and is due to the 
uncertainty on the flux calibration. The gap between -50 and 30 pc is due to our S/N cut 
which masked out the central, gas-deficient regions. 
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Fig. 4. — The surface brightness profile of NGC 1569 in the different narrow-band filters. 

Fluxes are corrected for Eg(B-V). The gap between -50 and 30 pc is due to our S/N cut 
which masked out the central, gas-deficient regions. 
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Fig. 5. — Panel a: the number of photons ionizing He [Q(He°)] is plotted as a function 
of the number of photons ionizing H [(^(H'')] as derived from the Br7 hne. The sohd hues 
trace the correlation between (5(H°) and (5(He°) at fixed total mass of the burst (10^, 10^ 
and 10^ Mq) and as a function of age (1, 3, 5, 7 and 10 Myr). The dashed hues show 
the correlation between (5(H°) and (5(He°) at fixed age and as a function of total mass 
of the burst. The black-filled circle represents the (5(H°) and (5(He°) values measured for 
NGC 1569 in this work. The errorbars are ilcr. Panel h: as before, but the solid fines 
here trace the predictions for a continuous star formation with different star formation rates 
(0.05, 0.10, 0.15 and 0.20 Mq yr~^). The line length corresponds to a time interval between 
1 Myr (bottom left corner of each fine) and 1 Gyr (top right corner); for ages older than 1 
Gyr (5(H°) and (5(He°) remain unchanged. 
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Fig. 6. — The star formation rate map of NGC 1569 derived from the dereddened Br7 image. 
The image size is 1.5 kpc x 1 kpc. The star formation rate is in units of M© yr~^ per pixel, 
and is displayed on logarithmic scale. 
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Fig. 7. — The azimuthal SFR density as a function of position along the semi-major axis, 
obtained from the dereddened Br7 image. The errorbar (± 0.11 dex) represents the system- 
atic uncertainty on SFR due to the uncertainty on the flux cahbration and color excess. The 
gap between -50 and 30 pc is due to our S/N cut which masked out the central, gas-deficient 
regions. 
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Fig. 8. — The 2D map of the [Fell] 1.64/im/Br7 flux ratio, after correction for reddening. 
The image size is 1.5 kpc x 1 kpc. 
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Fig. 9. — The azimuthal [Fell] 1.64/xm/Br7 flux ratio, corrected for reddening, as a function 
of position along the semi-major axis. The gap between -50 and 30 pc is due to our smoothing 
procedure which masked out the central, gas-deficient regions. 
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Table 1: Near-IR and optical total exposure times on source 



Filter 


Exposure Time (min) 


J 


20 


H 


55 


K 


121 


Hel 1.08/im 


30 


[Fell] 1.64/im 


67 


Br7 


111 


F606W 


326 


F658N 


77 



